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Abstract 

The incommensurate magnetic response of the normal-state cuprate perovskites is 
interpreted based on Mori's memory function approach and the t-J model of Cu-0 
planes. In agreement with experiment the calculated dispersion of the susceptibility 
maxima has the shape of two parabolas with upward and downward branches which 
converge at the antiferromagnetic wave vector. The maxima are located at (|, |±<5), 
(| ± 5, |) and at (| ± 5, ^ i 8), (j ± (5, 5 T <5) i n the lower and upper parabolas, 
respectively. The upper parabola reflects the dispersion of magnetic excitations of 
the localized Cu spins, while the lower parabola arises due to a dip in the spin- 
excitation damping at the antiferromagnetic wave vector. For moderate doping this 
dip stems from the weakness of the interaction between the spin excitations and 
holes near the hot spots. 
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One of the most interesting features of the inelastic neutron scattering in 
lanthanum cuprates is that for hole concentrations x > 0.04, low tempera- 
tures and small energy transfers the scattering is peaked at incommensurate 
momenta (|, ~ ± $), (~ ± 8, |) in the reciprocal lattice units 2ir/a with the 
lattice period a [1]. For x < 0.12 the incommensurability parameter 5 is ap- 
proximately equal to x and saturates for larger x [2]. The incommensurate 
response was observed both below and above T c [3]. Recently the analogous 
low-frequency incommensurability was observed in YBa2Cu307_ y [4]. This 
gives grounds to suppose that the incommensurability is a common feature 
of cuprate perovskites which does not depend on subtle details of the energy 
structure. However, for larger frequencies the susceptibility differs essentially 
in these two types of cuprates. In the underdoped YB^CusOr-,, both below 
and above T c a maximum was observed at frequencies u r = 25 — 40 meV 
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[5]. For this frequency the magnetic response is sharply peaked at the anti- 
ferromagnetic wave vector Q = (|, |). Contrastingly, no maximum at uj r was 
observed in lanthanum cuprates [6]. For even larger frequencies the magnetic 
response becomes again incommensurate in both types of cuprates with peaks 
located at (| ± 8, \ ± 5), (\ ± 5, \ =F 5) [5,7,8]. The dispersion of the max- 
ima in the susceptibility resembles two parabolas with upward and downward 
branches which converge at k = Q and uj uj r [8]. 

The nature of the magnetic incommensurability is the subject of active dis- 
cussion nowadays. The most frequently used approaches for its explanation 
are based on the picture of itinerant electrons with the susceptibility calcu- 
lated in the random phase approximation [9] and on the stripe domain picture 
[8,10]. In the former approach the low-frequency incommensurability is con- 
nected with the Fermi surface nesting which exists in the normal state or 
arises after the superconducting transition. In view of the similarity of the 
incommensurability in yttrium and lanthanum cuprates this approach implies 
the near resemblance of their Fermi surfaces and their change with doping 
which is highly improbable. Besides, the applicability of the picture of itiner- 
ant electrons for underdoped cuprates casts doubts. As for the second notion, it 
should be noted that in the elastic neutron scattering the charge-density wave 
connected with stripes is observed only in crystals with the low-temperature 
tetragonal or the low-temperature less-orthorhombic phases (La2- x Ba :r Cu04 
and L^-y-sNd^Sr^CuOzt) and does not observed in La 2 - 2; Sr :E Cu04 in the low- 
temperature orthorhombic phase [11]. At the same time the magnetic incom- 
mensurability is similar in these two types of crystals. It can be supposed that 
the magnetic incommensurability is the cause rather than the effect of stripes 
which are formed with an assistance of phonons. 

In the present work the general formula for the magnetic susceptibility derived 
in the memory function approach [12] is used. For the description of spin ex- 
citations in the doped antiferromagnet the t-J model of a Cu-0 plane is em- 
ployed. In this approach the mentioned peculiarities of the magnetic properties 
of cuprates are reproduced including the proper frequency and momentum lo- 
cation of the susceptibility maxima. The incommensurability for uj > uj r is 
connected with the dispersion of spin excitations. The incommensurability for 
lower frequencies is related to the dip in the spin-excitation damping at Q. For 
moderate doping this dip stems from the weakness of the interaction between 
the spin excitations and holes near the hot spots. Such form of the interaction 
constant follows from the fact that a decaying site spin excitation creates a 
fermion pair with components residing on the same and neighbor sites. 

In the memory function approach the imaginary part of the magnetic sus- 
ceptibility, which is directly connected with the cross-section of the magnetic 
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scattering, is described by the formula [13] 

/// k x = 4:fj? B u$sR(ku) 

X 1 J " [cu 2 - cu/ k 3?i?(M - cu 2 ] 2 + [cu/ k 3i?(kcu)] 2 ' 1 J 



where / k = (is k , -isi k ) _1 , (A, 5) = i / °° 5]), the angular brackets 

denote the statistical averaging, A(t) = e lHt Ae~ lHt with the Hamiltonian H, 
fi B is the Bohr magneton, s k = N^ 1 ^ 2 J2n e* kn s n with the number of sites 
iV and the z component of the spin s n on the lattice site n, is^_ = [s k , if], 
ujI = (is k ,-isi k )(4,si k )- 1 , R(kw) = -iJ™dte iuJt (A 2t ,Al) is the memory 
function with A 2 = « 2 s k — w k s k , 

^A 2t = (1 - P )(l - Pi)[A 2t ,H], A 2 , t=0 = A 2 , 
and the projection operators Pk, k — 0, 1, are defined as 

p fe g = (g,4)(A fc ,4)- 1 A fe 



with the operators A = s k and A\ = is k 



To describe the spin excitations of Cu-0 planes which determine the magnetic 
properties of cuprates the t-J model [14] is used. Its Hamiltonian reads 

H ^ ' ^nm"jnj"m(7 ~\~ 77 ^ ' ^nm^mi (2) 
nmcr nm 



where a n(7 = |ncr)(n0| is the hole annihilation operator, n and m label sites 
of the square lattice, a = ±1 is the spin projection, |n<r) and |n0) are site 
states corresponding to the absence and presence of a hole on the site, s n = 
\ J2a cr\na) (na| and s n = |ner)(n, — cr\ are the spin-i operators. Considering 
the interaction between nearest neighbor spins, described by the parameter J, 
and the hole hopping to nearest and next nearest neighbor sites, described by 
the parameters t and t', respectively, the terms in Eq. (1) are brought to the 
form [13] 



(isl, -is\) = 4(1 - 7 k)(J|C 1 | + tFx), 
col = lGJ^IdKl - 7k )(A + 1 + 7 k), 

^(M = ^x>kk</ du'Aiyj) 

k' — oo 

xA(k + V,co + .') nF{uJ + UJ,) - nF{uJ '\ 

UJ 
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where 7k = |[cos(A; x ) + cos(k y )], C\ = (s^s^jj and F 1 = (a^an+a) ar e the 
spin and hole correlations on neighbor sites, a is a vector connecting these sites, 
a ~ 1 is the vertex correction parameter, = 7k' — 7k+k' + 7 (7k' — 7k+k')> 
7 k = cos(k x ) cos(k y ), n F (u) = [exp(u;/T) + l] -1 with the temperature T, 
A(kcj) is the hole spectral function which is taken in the form 

AQuj) = ^ 2 (4) 



in this work. Here /x is the chemical potential, 77 is the artificial broadening, 
and £k is the hole dispersion. 

In Eq. (3), A oc £~ 2 where £ is the correlation length of the short-range 
antiferromagnetic order [15]. Thus, in the short-range order the frequency of 
spin excitations at Q is nonzero, in contrast to the classical antiferromagnetic 
magnons. As follows from Eq. (3), the dispersion of spin excitations has a local 
minimum at Q and can be approximated as c^k = [<^q + c 2 (k — Q) 2 ] 1 / 2 near 
this momentum [13]. 

To simplify calculations let us set T = 0. Notice that $sR(ku) is an even 
function of the frequency and consider the region u > 0. The integral in Eq. (3) 
reduces to J° w duo' A(k'uj') A(k + k', u + u/) and for the spectral function (4) 
is easily integrated. For r\ <C u the states with energies — uj < — fi < and 
< £k+k' — H < oj make the main contribution to this integral. 

In the following, we use the values of C±, F±, A and a calculated self-consis- 
tently in the t-J model on a 20x20 lattice for the range of hole concentrations 
< x < 0.16 [15]. The calculations were carried out for the parameters 
t = 0.5 eV and J = 0.1 eV corresponding to hole-doped cuprates [16]. For e k 
we apply the hole dispersion e k = -0.0879 + 0.55477 k -0.13277 k -0.01327 2k + 
0.0925 [cos(2k x ) cos(k y ) + cos(k x ) cos(2k y )] —0.02657^ proposed from the anal- 
ysis of photoemission data in Bi 2 Sr 2 CaCu208 [17]. Here the coefficients are in 
electronvolts. Analogous results can be obtained also with other model dis- 
persions [9,17]. 

The momentum dependence of x"(kcj) calculated with the above equations is 
shown in Fig. 1. As seen from this figure, there are three frequency regions with 
different shapes of the momentum dependence of x"(kc<j). The first region is the 
vicinity of the frequency ojq of the gap in the dispersion of spin excitations at 
Q. For the parameters of Fig. 1 ljq 37 meV. In this region x" is peaked at Q. 
For smaller and larger frequencies the magnetic response is incommensurate. 
The dispersion of maxima in x"(ku;) for scans along the edge and the diagonal 
of the Brillouin zone and their full widths at half maximum (FWHM) are 
shown in Fig. 2. 
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Fig. 1. Wave-vector scan of x" for T = 0, x « 0.12, /x = —40 meV, and ui = 70 meV, 
r/ = 3.5 meV (a), uj = 35 meV, rj = 3.5 meV (b), uj = 2 meV, 77 = 1.5 meV (c). 
Calculations were carried out in a 1200x1200 lattice. The solid lines correspond to 
scans along the edge of the Brillouin zone, k = (k, |); the dashed lines are for the 
zone diagonal, k = (k, k). 

To understand the above results notice that x" i n Eq. (1) has the resonance 
denominator which will dominate in the momentum dependence for uj > ujq 
if the spin excitations are not overdamped. Parameters of Fig. 1 correspond 
to this case. For uj > ujq the equation uj = uo^ determines the positions of the 
maxima in x"(kcj) which are somewhat shifted by the momentum dependence 
of ^sR (for simplicity the term with dtR is included in o; k ). Using the above 
approximation for uj^ we find that the maxima in x" are located near a circle 
centered at Q with the radius c~ 1 (uj 2 - ujq) 1/2 [13,18]. 

For uj < ujq the nature of incommensurability is completely different. It is 
most easily seen in the limit of small uj when Eq. (1) reduces to x"(kc<->) oc 
u;^ 4 9i?(kc<j) . u^ 4 is a decreasing function of the difference k — Q which acts 
in favor of a commensurate peak. However, if QR(kuj) has a pronounced dip 
at Q incommensurate peaks arise. For small x the dip appears due to the 
nesting of the ellipses forming the Fermi surface with the wave vector Q [19]. 
For larger x the mechanism of the dip formation is the following. For k = Q 
and small uj hole states which make the main contribution to ^sR are located 
near the hot spots (see the inset in Fig. 2). For these states the interaction 
constant qq^ in Eq. (3) is small which leads to the smallness of QR. With 
moving away from Q momenta of states contributing to the damping recede 
from the hot spots and ^k' an d QR grow. Thus, the damping has a dip at Q 
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Fig. 2. The dispersion of maxima in %"(ku;) for scans along the edge [k = (k, g), 
solid lines] and the diagonal [k = (k,k), dashed lines] of the Brillouin zone. The 
latter dispersion are drawn only in the frequency range in which these maxima are 
more intensive than those along the edge. Parameters are the same as in Fig. 1. 
Horizontal bars are FWHM for maxima along the edge. Inset: The Fermi surface 
for fj, = —40 meV (the solid line). Dashed lines show the boundary of the magnetic 
Brillouin zone, circles are the hot spots. 

which leads to the low-frequency incommensurability shown in Fig. lc. The 
smallness of gq k / near the hot spots is related to the fact that the decaying 
site spin excitation creates the fermion pair on the same and neighbor sites. 
As a result #kk' contains the functions 7k', 7k+k' and 7 k / — 7 k ' + k vanishing on 
the boundary of the magnetic Brillouin zone and at k = Q. 

The dispersion of peaks in x" which is similar to that shown in Fig. 2 was 
observed in YB^Cu^Oi^y and La 2 _ x Ba :E Cu04 [4,8]. As seen from Fig. 1, for 
frequencies u < uq the susceptibility is peaked at k = (|, |±<5) and (|±<5, §), 
while for lu > ujq the maxima are located at {\ ± 6, \ ± 6), {\ ± 5, \ T 6). This 
result is also in agreement with experimental observations [4,7,8]. Notice that 
in these crystals Fermi surfaces are similar to that shown in Fig. 2 [20]. The 
dependence of the incommensurability parameter S on x for small uo is shown 
in Fig. 3. In agreement with experiment (see the inset in Fig. 3) 5 grows 
nearly linearly with x up to x < 0.12 and then saturates. In these calculations 
the change of the hole dispersion with doping was not considered and the 
dependence 5(x) is connected solely with the variation of ujq oc oc x 1 / 2 
[15]. We found that the low-frequency incommensurability disappears when 
i] > uj. Besides, it disappears also if fi approaches the extended van Hove 
singularities at (0, tt), (tt, 0). In this case for k = Q the entire region of these 
singularities in which g k k' is nonzero contributes to QR. This may be the reason 
of the disappearance of the incommensurability in overdoped cuprates [2]. In 
the calculations the incommensurability disappears also in small lattices. 

As mentioned above, for the parameters chosen spin excitations are not over- 
damped near Q. As a consequence the frequency dependence of x"(Qc<j) has a 
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Fig. 3. The incommensurability parameter 6 vs. x for uj = 2 meV. The value of 5 for 
x = 0.043 was taken from Ref. [19]. Inset: experimental data [2] for L^-xSr^CuC^. 
Connecting lines are a guide to the eye. 



900 



600- 




co (meV) 

Fig. 4. The frequency dependence of \" . The solid lines are our results for T = 0, 
x » 0.12, n = -40 meV, rj = 3.5 meV, k = Q (a) and for k = (0.42,0.5) and 
the scaled hole dispersion (see text) (b). Squares are the odd susceptibility in the 
normal-state YBa 2 Cu 3 6 . 8 3 (x « 0.14 [21]) at T = 100 K and k = Q [5] (a) and 
the susceptibility in Lai 86Sro 14CUO4 for T = 35 K at the incommensurate peak [6] 
(b). 

pronounced maximum at uj m ujq which resembles the susceptibility observed 
in the superconducting and normal states of underdoped YBa2Cu 3 7 _j / . As 
seen in Fig. 4a, the experimental width of the maximum is approximately 
twice as large as the calculated one. Partly it is connected with the differ- 
ence in temperatures for the two sets of the data. Besides, the decrease of the 
hole bandwidth and the increase of the hole damping lead to a substantial 
growth of and to the overdamping of spin excitations [13]. An example 
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of such changes is shown in Fig. 4b where the results were obtained with the 
hole dispersion scaled by the factor 0.4. The calculations were carried out 
for k = (0.42, 0.5) which corresponds to the low-frequency incommensurate 
peak in Fig. lc. The growth of leads to the red shift of the maximum in 
x"(w). Its position is no longer connected with the frequency of spin excita- 
tions. The similar frequency dependence of x" is observed in L^-zSr^CuC^ 
[6] . The increased has no marked effect on the low- frequency incommensu- 
rability, however for the frequencies cjq > uj > 150 meV we found only broad 
commensurate peaks for the parameters of Fig. 4b. 

Our consideration was restricted to the normal state. Clearly the same mech- 
anisms lead to the incommensurate magnetic response also in the supercon- 
ducting state. 

In summary, it was shown that the calculated momentum and frequency de- 
pendencies of the imaginary part of the susceptibility x" i t ne dispersion and 
location of maxima in it and the concentration dependence of the incommen- 
surability parameter are similar to those observed in lanthanum and yttrium 
cuprates. The dispersion of the maxima in x" resembles two parabolas with 
upward and downward branches which converge at the antiferromagnetic wave 
vector Q and at the respective frequency of spin excitations. We relate the up- 
per parabola to the spin-excitation dispersion. The incommensurability con- 
nected with the lower parabola is related to the dip in the spin-excitation 
damping at Q. For moderate doping the dip arises due to the smallness of the 
interaction between spin excitations and holes near the hot spots. In agreement 
with experiment the incommensurate peaks which form the lower parabola are 
located at momenta (|, | ± 5), (| ± 5, |), while peaks in the upper parabola 
are at (| ± 5, \ ± 5), (\ ± 5, \ =F 5). Also in agreement with experiment the 
low-frequency incommensurability parameter 5 grows linearly with the hole 
concentration x for x < 0.12 and then saturates. By the variation of the 
hole bandwidth and damping the frequency dependencies of the susceptibil- 
ity which resemble those observed in YBa 2 Cu 3 7 _j / and La2_ x Sr x Cu04 were 
obtained. 

This work was supported by the ESF grant No. 5548. 
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